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Abstr act 

This paper describes Ihe process of analysis, design, digital 
implemenUtion and subsonic testing of an active controls 
flutter suppression system for a full span, frcc-lo roll wind 
tunnel model of an advanced fighter concept. The design 
technique employed a frequency domain representation of the 
plant and used optimization techniques to generate a robust 
multi-input/multi-output controller. During testing in a fixed- 
in-roll configuration, simultaneous suppression of both 
symmetric and antisymmetric flutter was successfully 
demonstrated. For a frcc-lo roll configuration, symmetric 
flutter was suppressed to the limit of the tunnel lest envelope. 
During aggressive rolling maneuvers above the open -loop 
flutter boundary, simultaneous flutter suppression and 
maneuver load control were demonstrated. Finally, the flutter 
suppression controller was reoptimized overnight during the 
test using combined experimental and analytical frequency 
domain data, resulting in improved stability robustness. 

Introduction 

An advanced fighter aircraft design which exploits, rather 
than avoids, wing flexibility to provide improved aerodynamic 
performance is likely to require an active flutter suppression 
system (FSS) to remove dynamic structural instabilities 
(particularly if FSS is employed in the design process to 
minimize weight). If the FSS is required for stabilization 
within the operational envelope, it is essential that proper FSS 
functioning be maintained during aggressive maneuvers. This 
paper describes the design and wind-tunnel lest of an active 
FSS for a configuration that exploits wing flexibility. 
Reliability aspects arc not addressed. 

The test vehicle used in the study was the Active Flexible 
Wing (AFW) wind-tunnel model built by Rockwell. 1 Testing 
was conducted in the NASA Langley Research Center 
Transonic Dynamics Tunnel (TDT). The model (sec fig. 1) 
was sting mounted and could be rolled about the sting axis 
(frce-to-roll configuration) between plus and minus 
145 degrees. Roll motion could be prevented through the use 
of a mechanical pin prior to testing, or stopped during a 
maneuver through the use of a hydraulically actuated roll 
brake (fixed- in -roll configuration). 

Active controls flutter suppression of the AFW wind-tunnel 
model was tested during TDT entries, in 1989 and 199 1, 2 
using a dedicated programmable digital controller. 3 For the 
1989 entry, only the fixed-in-roll configuration was tested for 
plant identification and for flutter suppression. 4 * 5 Results in 
the present paper pertain to the 1991 entry, the primary 
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Fig. 1 Sketch of AFW wind-tunnel model. 


objective of which was to demonstrate simultaneous 
application of digitally implemented multi-input/multi-output 
(MI MO) flutter suppression and maneuver load controllers 
while performing aggressive rolling maneuvers above the 
open -loop flutter boundary. 

Flutter suppression testing was conducted subsonically in 
air, operating at atmospheric static pressure. For the fixed-in- 
roll configuration, both symmetric and antisymmetric flutter 
were predicted to occur within the TDT operating capability. 
For the frce-to-roll configuration, only symmetric flutter 
(unaffected by roll freedom status) was predicted to occur 
within the TDT operating capability. Thus, one control law 
was developed to suppress symmetric flutter and another was 
developed to suppress fixcd-in-roll antisymmetric flutter. Both 
control laws were active, operating in parallel, for fixed-in-roll 
flutter suppression testing. Only the symmetric law was active 
for the frce-to-roll tests. 

Modeling 

Linear mathematical models were generated for design and 
preliminary evaluation of candidate controllers. The 
Interaction of Structures, Aerodynamics and Controls (ISAC) 
system of programs 6 * 7 * 8 played a key role in providing linear 
models (see fig. 2). ISAC received inputs that included in 
vacuo modal data, planform geometry, aerodynamic paneling 
specifications, turbulence power spectral density 
characterizations, sensor locations, and actuator transfer 
function descriptions. It then provided linear aeroelastic 
equations of motion at specified conditions for use in analysis 
and design. A frequency domain form of the equations of 
motion made direct use of tabular unsteady aerodynamic 
forces computed using a doublet lattice code contained within 
ISAC. A finite dimensional stale space form employed 
rational function approximations of the unsteady aerodynamic 
data. 8 
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Fig. 2 Modeling via ISAC. 

Closed-loop operation of the candidate controllers was 
validated prior to wind-tunnel testing with simulation analyses. 
ISAC provided intermediate, linear, mathematically generated 
data necessary for the creation of the simulations (see fig. 2). 
The simulations then incorporated separate left and right 
models for the actuators, actuator deflection and rale limits, 
and quantization effects for the analog/digital conversion 
process. 9,10,11 Stale space models could be extracted from 
the simulation for control law design purposes by perturbing 
the simulation model. The simulation allowed both 
symmetries to be run simultaneously in the presence of 
independently generated symmetric and antisymmetric 
turbulence. 

The frequency domain form of the equations of motion was 
employed in development of the controller design discussed 
herein. This form of the equations allowed numerical 
compulation of a frequency response for any output/input pair 
and direct replacement with experimental data when they 
became available. The frequoncy domain form of the 
equations of motion was presented in reference 4. 

Modal data from a structural model developed prior to the 
1989 TDT entry were employed in generating the equations of 
motion used for analysis and design. This selection was made 
in preference to a post- 1989 model that was also available 
because the character of the resulting analytical frequency 
response predictions was in closer agreement with experiment, 
and both models exhibited comparable levels of inaccuracy in 
predicting the flutter frequencies. 


Flutter and Inpui/Qutput Characteristics 

Analyses were made, using the linear mathematical models, 
of the predicted stability and response characteristics of the 
AFW wind-tunnel model in order to determine the flutter 
characteristics and to assess candidate controls and sensors to 
use in the control laws. The trailing edge outboard (TEO) and 
trailing edge inboard (TE1) pairs of control surfaces (see fig. 1) 
were effective for flutter suppression. Each pair of wing 
accelerometers was characterized in terms of its location, the 
leading edge outboard (LEO), trailing edge inboard (TEI), 
trailing edge outboard (TEO), and wingtip (TIP) pairs. Key 
analysis results are presented below. 

Figure 3 presents loci of symmetric poles generated by 
varying dynamic pressure, expressed in pounds per square foot 
(psf), from 0 psf to 325 psf in 25 psf increments. The analysis 
was performed using a linear state space representation that 
contained single-lag rational function approximations of the 
unsteady aerodynamic forces. The loci of poles depicts a 
primary flutter mechanism that is an example of classical 
airplane flutter, with approximate coalescence (near 70 radians 
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Fig. 4 Frequency responses, 
(symmetric, open -loop, 175 psf) 



Fig. 3 Root locus. 

(symmetric, open-loop, dynamic pressure variation) 


per second) of the damped frequencies of two roots 
corresponding to modes that, in vacuo, were predominantly 
first wing bending and first wing torsion. The predicted flutter 
dynamic pressure was 248 psf. A second pole, near 
230 radians per second, was approaching instability at 325 psf, 
indicating the presence of a secondary flutter mechanism. The 
significant, although smaller than desirable, frequency 
separation between the primary and secondary flutter 
mechanisms facilitated the task of primary flutter suppression 
without secondary flutter degradation. 

Figure 4 presents information typical of that used to select 
control surfaces and sensors to employ in the controller. Bode 
plots are presented for two of the symmetric components 
( ^hf/SteOc and 2 leo/5tei c ) of the AFW element shown in 
figure 5, at a dynamic pressure of 175 psf. Inspection of this 
figure indicated that Eleo/^TEIc had die undesirable properly 
of responding to modes in the higher, secondary flutter 
frequency range. The undesirable high frequency response 
was present for all sensors other than those at the wingtip. The 
£rip/§TEO c transfer function exhibited desirable high 
frequency rolloff. However, the primary flutter mode was 
more in evidence for 2leo/&TEI c than for 2 tip/&teo c - ^h e 
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relative prominence of the flutter mode among the various 
transfer functions was strongly influenced by the variation in 
location of a critical, lowly damped zero relative to the pole 
that became unstable (see reference 1 2 for a discussion of the 
impact of the critical zero upon controller design). The 
analysis indicated that all of the wing accelerometers, with 
appropriate high frequency filtering, were viable candidates for 
feedback. 

For the fixed -in-roll configuration, the loci of antisymmetric 
poles (not shown) exhibited behavior similar to that of figure 3 
with predicted flutter at about the same frequency and with 
onset («252 psf) that was higher in dynamic pressure by only 
about 1 .6 percent. Likewise, the frequency separation 
between Lhe primary flutter mechanism and a cluster of higher 
frequency modes was similar to that of the symmetric case. 

Pretest Design Activity 

The design philosophy adopted was to seek an unscheduled, 
MIMO flutter suppression controller. The design was 
accomplished in the continuous domain although the controller 
was implemented digitally using a Tustin transformation. A 
number of design considerations were addressed. 

Design. Considerations 

Three dynamic elements that will recur in the discussion of 
controller dynamics (a general second -order element, a first- 
order highpass clement, and a first-order lowpass clement) are 
defined here as 

e = e(^n, C<1. coj) 

= (s 2 +2 ; n ®n S + (On 2 ) / (s 2 +2 G) d S + <0 d 2 ). (1) 

HP = HP(a) = s / (s + a), (2) 

LP = LP(c) = c / (s + c), (3) 

where s is the Laplace variable, £ n and arc numerator and 
denominator damping ratios, respectively, co^ and co<] arc 
numerator and denominator natural frequencies, respectively, 
a is the highpass filter pole, and c is the lowpass filter pole. 

The frequency separation between the roll dcgrcc-of- 
freedom and the flutter mechanism was sufficiently large to 
conclude that FSS and load-limiting rolling maneuver control 
laws could safely be designed separately and then operated 
simultaneously provided filtering was included to prevent 
interference. Accordingly, highpass dynamic elements (HP, 
eq. (2)) were included in the FSS controller to ensure low 
frequency separation by attenuating any low frequency 
disturbance to or from the load -limiting active control law. A 
lowpass filter (LP, eq. (3)), with break frequency well below 
10 Hz, in the maneuver load control (MLC) system would 
complete the frequency separation objectives. 

Figure 5, which refers to either of the symmetries of motion 
(symmetries were decoupled in the linear models), indicates 
that all four accelerometers on each wing were used in the 
feedback control law to drive both trailing edge control 
surfaces. Multiple sensors and control surfaces were used both 
to satisfy program MIMO objectives and to exploit the 
favorable characteristics of the various sensors and controls 
that were seen in the analysis discussed previously. 
Attenuation of undesirable high frequency sensor signal 
components in accelerometers other than £pip was 
accomplished using available analog notch filters, Np(s) = 
c(0.08, 0.32, 2 k (32), 2k (32)), sec eq. (1). 

The sample rale for digitally implemented controllers was 
prescribed to be 200 Hz. To reduce aliasing, a first-order 
analog filter with break frequency of 25 Hz (157 rad/scc) was 
incorporated on each channel to be sampled. These 
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Fig. 5 Design plant definition. 


antialiasing fillers also provided some attenuation between 25 
and 100 Hz and, thereby, reduced signal strength in the range 
of higher frequency structural modes. Phase lag at 10 Hz (near 
the flutter frequency) due to this filter was about 22 degrees. 

The plant used in design of the continuous controller (see 
fig. 5) contained analog notch and antialiasing filters and an 
additional time delay element that was included to compensate 
for effects of digital implementation. Since, in actual digital 
implementation, some time would be required to generate 
control command outputs after receipt of sensor inputs, the 
control laws were implemented in the digital controller such 
that control commands generated based upon sensor inputs 
received at one sampling instant were held until the next 
sampling instant. Therefore, a net delay of 1-time step 
associated with holding the controller output was included in 
the "design plant” model. This amounted to a phase lag of 
1 8 degrees at 1 0 Hz. 

Use of an unscheduled controller required closed-loop 
stability over the full test range of dynamic pressures (100 psf- 
to-325 psf) based upon the analytical simulation model. Gain 
margins of ± 4 dB and phase margins of ±30 degrees, or their 
multivariable equivalents, were also required over this range. 
In addition, stability in the presence of ±10 percent variations 
in the frequencies of the coalescing modes was required at a 
300 psf evaluation point. A design point of 325 psf was 
chosen. 

Actuator rale saturation can effectively induce lag and 
reduce the amplitude of control surface deflections. At 
325 psf, the open-loop time-to-double amplitude was predicted 
to be about 1/10 of a second. For this level of instability, 
actuator rale saturation of a pair of actuators for even a brief 
period of time in response to wind-tunnel turbulence could 
cause unacceptably large growth of the flutter mode. This 
reinforced the need for restriction on actuator rate 
requirements. The TEO surfaces were conservatively assumed 
lo have a no-acrodynamic load rate limit of 150 degrees per 
second. It was assumed that no rate saturation for a 3-standard 
deviation turbulence velocity magnitude was adequate for 
assuring that actuator rate saturation was sufficiently unlikely. 
This constrained the RMS rate for a 1 -standard deviation 
lurbulcncc velocity to be less than 50 degrees per second. 
Definition of the assumed turbulence characteristics is 
presented in reference 1 1 . 

Design Approach 

In this section the approach taken to determine the 
controller, K(s), is outlined. The approach was the same for 
each symmetry and separate designs were obtained. The first 
step was lo explicitly specify the controller structure, including 
all dynamic elements, thereby displaying selectable design 
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variables. Optimization techniques were employed to select 
values for the design variables. The optimization process look 
place in Iwo distinct stages. For each stage, a conjugate 
gradient algorithm 13,14,15 W as used to minimize a cost 
function representing dissatisfaction with respect to a set of 
design criteria. For the first stage a dominant control 
surface/sensor pair was selected. The trailing edge outboard 
control surfaces and the accelerometer pair nearest the wing 
tips were selected because the combination exhibited desirable 
high frequency rolloff characteristics and adequate response in 
the flutter frequency range at the 325 psf design point. 
Parameters determining the controller dynamics were then 
optimized for this single-input/single output (SISO) 
compensator to increase tolerance to multiplicative and 
additive plant errors and to increase tolerance to errors in the 
predicted flutter frequency. The second stage, employing 
fixed dynamics from the first stage, determined the 
coefficients used for blending the four pairs of sensors into one 
composite sensor and for distributing the filtered feedback 
signal to the two pairs of controls. The choice of blending and 
distribution coefficients was made so as to achieve MI MO 
robustness to additive plant errors and to errors at plant input 
and output. 

SISO Optimization 

The form of the SISO compensator for the first stage 
optimization is shown in figure 6. There was an overall gain, 
(kl), a highpass element (HPj(ai)) that has already been 
discussed, a broad notch (Np = e(£ n N, £djsi» W N. ®n)) that 
reduced the observability of the coalescing pole that became 
more damped, a bandpass (amplification) clement (Aj; = 
c (CnA* <*>a)) partially removed the negative 

impact of a critical zero in the flutter frequency range (refs. 4 
and 12), a higher frequency notch to avoid adversely affecting 
higher frequency modes, and a lowpass element for rolloff. 
The higher frequency notch and lowpass elements were unity 
at this stage and will be discussed below, where first used. 

SISO Design and Cost Function Variables, The eight 
design variables used in the SISO optimization were 

0^1* a t» ’CdN* CriA* ®a)- 
The highpass frequency was constrained to be between 10 and 
25 rad/sec, and all damping ratios were constrained to be 
between 0.1 and 1. These design variable constraints were 
enforced through a trigonometric mapping technique. 16 

The dependent variables used to form the cost function 
were: 

&o = min { abs[ 1 +f 0 (ico) ] } 

0) 

6p = min [ abs[ 1 + f 0 (ito) ] / abs[ k(ico) ] }, g/deg 
to 

I = f True ^ closed loop stable 

L “ i False if closed -loop unstable 

where go varied from 10-to-400 rad/sec (1.6-to 64 Hz). The 

symbols p(ito), k(ico), and f 0 (ico) = p(ico) k(ico) = k(ico) p(ico) 

represented frequency responses for the SISO plant, controller 
and loop-transfer functions, respectively (sec fig. 6). The 
variables o 0 (a>) and a p (co) were singular values which 
represented tolerance to unstructured multiplicative error and 
to unstructured additive plant error, respectively. The ( A ) 
represented the minimum over the range of frequencies 
considered. The logical variable, L, indicated whether the 
closed -loop system was stable. Closed -loop stability was 
computed by knowing the stability of the open loop plant and 
requiring the appropriate number of counterclockwise 


encirclements of the critical point (-1) for a polar (complex 
plane) plot of F 0 (ico). 

Terms were included in the cost function to desensitize the 
controller to errors in the predicted frequencies of the poles 
that coalesced to produce flutter. Design criteria were 
specified not only for the nominal plant, but also for eight 
variational plants obtained by perturbing the frequency of at 
least one of the two coalescing poles by ±10 percent. The 
variational plants were not generated. Instead, an alternate, 
approximately equivalent approach was adopted: namely 
perturbing the frequencies of related compensator variables. 
The frequency of the amplification element was perturbed to 
mimic a corresponding frequency shift in the opposite 
direction for the plant pole that was unstable. Likewise, the 
frequencies for the notch and the highpass elements were 
perturbed in unison to mimic a shift in the frequency of the 
stable coalescing plant pole. Each iteration of the optimization 
required evaluation of the cost function components for the 
nominal compensator and its eight variations. 

During early stages of the study, RMS control surface rates 
had been computed for inclusion in the cost function, but were 
removed from the optimization after observing that they rarely 
exceeded the 50 dcg/scc design requirement. 



a 1 co N co A co 1 c 
co (rad/sec) 


Fig, 6 SISO controller dynamics. 

SISO Cost Function and Component Scaling . The design 
variables were scaled to improve the mathematical 
conditioning of the optimization, but proper scaling of the 
dependent variables used in forming the cost function was 
more important. The scaling of these cost function 
components had a direct impact on the solution achieved. The 
cost function had the form: 

9 

Jsiso = Vo , I * * 4 + V p] 4 + E { Vo ,,, 2 + VjJ } (4) 

m=2 

where m was an index indicating either nominal (m = 1) or 
perturbed (m = 2-to-9) compensator dynamics. V 0 and V p 
were sealed violations of constraints, defined as follows: 

If was true (stability) 

Vo™ = S 0 max(0, T^-d^), 

V Pm = S P maxIO.TpH, -d^), 
otherwise 

V om - U o(To, n + A 0m ). 
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Vp„, = Up (T^ +&PJ, 

where 

S 0 = 10, Sp = 1 0 deg/g were scale factors if stable, 

U 0 = 14, Up = 14 deg/g were scale factors if unstable, 

T 0 | = 0.9, Tpj = 0.9 g/deg were target min. sing, values, 

T om = 0.6, T ^ = 0.6 g/deg were target min. sing, values, 
for m > 1 . 

Some comparisons can be made between nominal (m = 1) and 
variational (m = 2-io-9) terms in the cost function (sec eq. (4)). 
The nominal terms were emphasized by selection of larger 
target minimum singular value magnitudes (0.9 vs 0.6). The 
selection of unachieved 0.9 magnitudes caused the nominal 
scaled violations V ol and V pl to exceed 1 throughout the 
optimization. Consequently, raising these nominal scaled 
violations to the fourth power in the cost function, as 
compared with only squaring the variational sealed violations, 
further accentuated the nominal case. 

Equal target minimum singular value magnitudes and equal 
scale factor magnitudes were chosen for plant additive error 
and loop transfer function multiplicative error despite the 
dimensional incompatibility. This choice was made because, 
for prototype stabilizing SISO compensators, the smallest 
values for a p (to), in g/deg, and for a 0 (tu) were of comparable 
size. Different target magnitudes and different scale factor 


Table 1 Tolerance to flutter frequency and gain variations 
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magnitudes would generally be required for variables having 
disparate magnitude ranges. 

Cost function component definitions changed at the stability 
boundary such that singular values were driven smaller rather 
than larger for the closed-loop unstable condition. Thus, the 
requirement for an initially stabilizing selection of parameters 
for the nominal controller and all of its variations was relaxed. 
Furthermore, scale factor changes produced a discrete jump in 
the cost at the stability boundary. A scale factor discontinuity 
was acceptable in this case because the stability boundary was 
not a satisfactory solution. As a result of the discrete jump and 
the composition of the performance function, once a 
stabilizing controller was achieved, destabilizing controllers 
were removed from consideration. 

SISQ- Performance Predictions . Table 1 summarizes the 
results of the SISO optimization for the symmetric control law 
in terms of closed-loop insensitivity to frequency variations in 
the open-loop plant. The closed-loop system was stable at 
300 psf (the required evaluation point) for all variational 
combinations of the frequencies, and had 4 dB or more gain 
margin for all combinations except one, as indicated. 

Figure 7 presents Nyquist plot evaluations of the 
antisymmetric SISO nominal and variational controllers at 
300 psf (open-loop unstable, critical point at -1). The solid 
curve is for the nominal case. This figure shows that all 
variational controllers arc stable with gain margins exceeding 
± 4 dB. The figure also contains the information needed for 



Real 

Fig. 7 Nyquist plots for SISO controllers, 
(fixcd-in-roll, antisymmetric, 300 psf) 

determination of phase margins, minimum singular value and 
maximum loop gain for each of the nine controllers. 

M1MO Optimization 

The second stage of the optimization pertained to the 
MIMO controller. As shown in figure 8, four pairs of sensors 
and two pairs of control surfaces were employed. The 
dynamic elements Ni, N 2 , and LP were unity at this stage. 
They will be discussed below, where first used. The dynamic 
elements found during the SISO optimization were retained, 
unchanged, in the $teo c / ztip channel. However, for the 
channel containing information from the other three sensor 
pairs, the highpass break frequency (a2 in HP 2 ) was increased 
as compared to sl\. The selection of a 2 was made such that, at 
the flutter frequency, the additional lead compensated for the 
lag due to the 32 Hz analog notch filter (see Figure 5) on each 
of the three pairs of inboard sensor channels. This produced 
phase matching at the flutter frequency for the primary and the 
notched sensor channels, at the expense of reduced phase 
margins for the notched channels. 

MIMO Design and Cost Function Variables . The design 
parameters for the MIMO optimization were the overall gain, 
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Fig. 8 MIMO controller structure. 







(Jez), the blending coefficients For the inboard sensor pairs (bi t 
b2» bj), and the distribution coefficient (d) for the inboard 
control surface pair. These five parameters are shaded in 
figure 8, No constraints were imposed on the design 
parameters. 

The dependent variables used to form the cost function 
components were: 

& s min {a f 14x4 + P(ico) K(io>) 1 } 

G) 


&C S min {a f 12x2 + K(ico) P(iw) ] } 

<D 


&p 3 min { 1 /a [ K(ico) [ 14x4 + P(ico) K(ico) ] 1 ] } , g/deg 

CO 



0 if closed-loop stable 
100 if closed-loop unstable 


where af ] denoted the operation of finding the singular values 
of the matrix within the brackets, (_) and (“) denoted 
minimum and maximum, respectively, and, as in the SISO 
case, ( A ) denoted the minimum over <0. The indicated 
minimization over (o covered the frequency range from 10- to- 
400 rad/sec (1.6-to-64 Hz). Symbols P(ico) and K(ico) 
represented frequency response matrices for the 2- input, 
4-output plant and for the 4-input, 2-outpul compensator, 
respectively (see fig. 8). The functions £s(o>) and Qc(u>) were 
minimum singular values of plant output and input return 
difference matrices which represented tolerance to 
unstructured multiplicative error at the sensors and at the 
controls, respectively. 17 The function 2p(to), in g/deg, was the 
reciprocal of a maximum singular value and represented 
tolerance to unstructured additive plant error. 

The variable B was a penalty to be added to the cost 
function if a candidate set of design variables caused closed- 
loop instability. Stability was checked at each iteration after 
evaluating the frequency response, f ss (ico), of the loop transfer 
function at the single strand point identified in figure 8. 
Closed-loop stability was computed by knowing the stability 
of the open-loop plant and requiring the appropriate number of 
counterclockwise encirclements of the origin for a polar plot of 
[l+f ss (io>)]. (Although this was a reliable indicator of absolute 
stability, stability margins defined at the single strand point 
were potentially nonconservalive 4 and, therefore, were not 
used). 

MTMO Cost Function and Component Scaling . Values of 
the design variables were sought which minimized the 
following cost function (only fixed, nominal compensator 
dynamics were considered) 

Jmimo = V s 2 + V t 2 +V p 2 + B, 
where 


V s = S s max(0,T s -&j) 

V c ” Sc max( 0, T c -&c) 

V p = S p max(0,T p -&p) 

where 

S s = S c = 10, and S p = 10 deg/g were scale factors, and 
T g = T c = 0.6, and T p = 0.6 g/deg were target min. sing.vals. 
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As in the SISO case, selection of equal scale factor magnitudes 
and equal target lower bound magnitudes for each of the 
dimensional and nondimcnsional variables was based upon the 
various minimum singular values having approximately the 
same magnitudes. 

Since the MIMO optimization used zero starting values for 
the blending coefficients for the additional sensors and 
controls, and since the initial SISO controller was known to 
stabilize the closed-loop system, no provision was made to 
drive singular values smaller for a closed-loop unstable 
situation. Rather, consideration of closed-loop unstable 
controller coefficients was precluded by making the stability 
boundary approach the characteristics of a hard constraint 
through the use of the large discrete jump (B) in the cost 
function. 

For each symmetry, the controller that was obtained as a 
result of the MIMO optimization process was predicted to 
suppress flutter over the anticipated test range of dynamic 
pressures, to satisfy all robustness criteria, and to call for RMS 
actuator rales that were well within the 50 deg/sec design 
requirements. 

MIMO Simulation and Controller Adjustment . Candidate 
control laws designed based upon linear models at 325 psf 
were evaluated over the anticipated test range of dynamic 
pressures using the nonlinear simulations. The closed-loop 
system was stable at all points for both fixed and free-to-roll 
configurations. However, the gain (k2) was adjusted more 
toward the center of its stable range. In addition, notch filters 
(Ni) were added at the single strand point (see fig. 8) to reduce 
transient, turbulence induced responses observed in the 20 Hz- 
to-50 Hz frequency range. The notches added (which were to 
be implemented digitally after the prewarping discussed 
below) were Nt(s) = e(0.40, 0.60, 2 n (23.9), 2 n (23.9)) 
symmetrically and Ni(s) = c(0.20, 0.60, 2 n (50.0), 2 n (50.0)) 
antisymmetrically. New experimentally determined actuator 
transfer function estimates were incorporated which, 
fortuitously, approximately compensated for the phase lag due 
to the added notches, removing the need for reoptimization. 

Figure 9 shows time histories generated using the 
simulation from reference 9 for a candidate pretest control law 
that did not have the notches, Nj. The results represent a 
flxcd-in-roll condition, with no turbulence input, at the design 
point of 325 psf. Beginning at 0.1 second, a 10 Hz 
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Fig. 9 Nonlinear simulation performance, 
(fixed -in -roll, 325 psf, no turbulence) 




Table 2 Parameter values for final controller 



Symmetric 

Antisymmetric 

HPt 

ai = 10.0 1/scc 

a] = 10.0 1/sec 

Nf 

Cn N = 0-20 

;„ N = 0.14 


Cd N = 0.57 

Cd N = 0.33 


C0n N = 47.1 rad/scc 
(OdM = 44.6 rad/scc 

con = 44.0 rad/sec 

Af 

Sn A = 0.55 

Cn A = 0.37 


Cd A = 0.12 

W = 0.12 


con A = 68.5 rad/scc 
a>d A = 64.7 rad/scc 

«a = 69.7 rad/scc 

Ni a 

= 0.40 

Cni = 0.20 


^d, = 0.60 

£d, = 0.60 


®l = 150 rad/scc 

(Ol = 314 rad/scc 

LP 8 

c = 364 1/scc 

not used 

hp 2 

a^ = 23.2 1/sec 

a 2 =18.2 1/scc 

n 2 * 

Cn 2 = 0.16 

not used 


= 0.48 

not used 


(02 = 251 rad/scc 

not used 

Gain Coefs 

k2 = 0.697 deg/g 

k 2 = 0.675 deg/g 


b[ = 0.586 

b[ = 0.331 


b 2 = 0.006 

b 2 = 0.213 


b 3 = 0.265 

b 3 = 0.437 


d = 0.453 

d = 0.552 


a prior to prcwarping. 


antisymmetric doublet command of 1 degree magnitude was 
input to the TEO actuators. The antisymmetric component of 
the resulting (closed-loop) control surface motion, 5rFO aS y. is 
shown in the figure. The antisymmetric wingtip accelerometer 
response, 2rip asy , was damped out in approximately 2 cycles 
after the command ended at 0.2 second. Small effects due to 
actuator asymmetries were apparent in the symmetric wingtip 
accelerometer response, £np sym . The small persistent 
oscillations in all responses were ihe result of a limit cycle 
arising from quantization effects. 

Controller Implementation 

Digital Implementation of Continuous Design 
The control laws (all designed in the continuous domain), 
were implemented digitally, both in the simulation and in the 
digital controller, using a Tustin transformation with a 200 Hz 
sample rate. Prcwarping was applied for those dynamic 
elements which required an analog/digilal frequency response 
match at a frequency above 30 Hz. The prcwarping was 
nonstandard. An iterative procedure was employed wherein a 
"desired” analog transfer function was specified, and a 
"modified" analog transfer function of the same form and order 
was sought which had a digital transform with frequency 
response close to that desired at and below a chosen frequency 
for a match. 

Frequency Shift in Test Controller 

A frequency shift was made in key controller dynamic 
elements for implementation during the lest. The shift was 
based upon data gathered during the 1989 lest which showed 
that the 1989 analytical model of the plant ovcrprediclcd the 
frequencies for the coalescing modes at flutter, for both 


symmetries, by nearly 2 Hz. It had also been shown that the 
predicted phase at the predicted frequency of peak response 
was in close agreement with the observed phase at the 
observed frequency of peak response. 4 The amplitude of 
accelerometer response to control deflections was also 
somewhat ovcrprediclcd. Table 1 indicates that, for the 
symmetric SISO controller, the situation where the actual 
frequencies and gain were lower than predicted was the least 
well tolerated combination of errors. For these reasons, the 
notch clement (Np) and the amplifier element (Ap) of the 
controller were shifted downward by 2 Hz for both symmetries 
prior to implementation for test. 

Testing and Controller Updates 

Impact of Early Testing on FSS 

Early in the entry, investigators encountered an unexpected 
64 Hz buzz of the trailing edge outboard control surfaces when 
the FSS loops were closed. New servovalves were installed 
for each actuator and an analog 64 Hz notch was inserted in 
each trailing edge outboard command channel. The buzz was 
removed, but approximately 13 degrees of phase lag at 10 Hz 
resulted. 

Plant estimates were obtained below flutter during 
subcritical open-loop tests, and above flutter by analysis of 
results from post critical closed-loop tests of other FSS designs 
that had been tested earlier in the entry. The frequency -shifted 
controller designed prior to test was predicted, based upon the 
plant estimates, to be stabilizing. However, phase adjustments 
were made prior to closed -loop tests to improve robustness 
characteristics. The phase adjustments were made by 
modifying the highpass break frequencies and incorporating a 
lowpass filter, LP(c) (see figure 8, eq. (3) and table 2). The 
resulting symmetric law was then successfully tested to the 
tunnel limit in the free-to-roll configuration, reaching a 
dynamic pressure 25 percent above the observed (»235psf), 
open-loop flutter dynamic pressure. 

Overnight Controller Reoptimization 

An overnight reoptimization of the symmetric controller 
was accomplished during the wind-tunnel test using a hybrid 
of experimentally derived frequency responses for frequencies 
between 5 and 20 Hertz, and analytically predicted frequency 
responses for frequencies between 20 and 64 Hertz, all for a 
285 psf condition. The frequency content for the excitation 
signal for the control surfaces used for feedback was limited to 
20 Hertz and below, so that experimentally derived frequency 
responses above 20 Hertz were not available. An additional 
40 Hz notch (N2) was implemented digitally (after 
prewarping) on the output of HP2 pertaining to the 
accelerometers other than the tip accelerometers (see figure 8 
and table 2) to further safeguard against potential aggravation 
of the stability and response of higher frequency modes (the 
reoptimized controller utilized these inboard accelerometers 
more than the symmetric design based upon the analytic model 
did). 

The SISO portion of the reoptimization did not include the 
variational terms (m=2-to-9 in eq. (4) ). Two additional design 
variables were also included by allowing distinct numerator 
and denominator natural frequencies in the notch (Nf) and 
amplifier (Ap) dynamic elements. As mentioned previously, a 
2 Hz frequency shift had been made in the controller dynamics 
in the flutter frequency region based upon data from a previous 
lest. The reoptimization made only small additional shifts in 
the controller frequencies that had been selected prior to test. 
The parameters for the reoptimized controller, including 
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Frequency (Hz) 


Fig. 10 Minimum singular values for reoptimized control law. 
(symmetric, 285 psf) 


configuration. A maneuver load control (MLC) roll rate 
tracking system 18 was also active but was commanding zero 
roll rate. Figures 11 and 12 show absolute stability and 
stability robustness assessment information obtained in near- 
real-time during the test through the use of a Controller 
Performance Evaluation (CPE) analysis procedure. 5 - 19 These 
results were obtained using MIMO extensions to Fast Fourier 
Transform based S1SO controller performance evaluation 
techniques. 20 * 21 

Figure 1 1 shows MIMO Nyquist plots which, when 
accompanied by knowledge of open-loop stability 
characteristics, provide a definitive assessment of closed-loop 


definition of all notches to be implemented digitally (after 
prewarping) are presented in table 2. 

The reoptimization of this symmetric control law was 
accomplished overnight. In-flight excitation of the control 
surfaces and use of a computer that was two orders of 
magnitude faster would have allowed the plant frequency 
response estimates and controller redefinition to be completed 
quickly enough to perform a tuning type adaptation of the 
controller to changes in plant characteristics arising from slow 
(relative to the adaptation process) Mach number and dynamic 
pressure variations. 

Figure 10 shows the predicted reoptimized controller 
performance based upon the hybrid experimental and 
analytical data. The local minima seen in the various 
minimum singular value curves arc relatively wel] balanced in 
the O-to-20 Hz frequency range with the smallest values of 
approximately 0.4 representing a significant level of 
robustness. 

Test Results with Reoptimised Controller 

The reoptimized controller was successfully tested with the 
AFW wind-tunnel model in both its frec-to-roll and fixed-in- 
roll configurations. Results will now be shown that provide a 
more detailed description of the performance of the 
reoptimized controller. 

With the reoptimized controller employed, symmetric flutter 
was again suppressed to the tunnel limit in the frce-to-roll 






Frequency (Hz) 

c) dynamic pressure of 290 psf. 


Fig. 1 1 Experimental MIMO Nyquist plots, 
(symmetric, reoptimized FSS with MLC) 
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Fig. 1 2 Experimental singular values, 
(symmetric, reoptimized FSS with MLC) 







stability. The MIMO Nyquist plots are polar plots of 
del[ 12x2 + K(ia») P(ito) ] = del[ 14 x 4 + P(ico) K(ico) ] 
or, equivalently, for a rank one loop transfer matrix, 

= 1 + f M (io>). 

The open -loop system was stable at the two lower dynamic 
pressures and unstable at the two higher dynamic pressures. 
Thus, figure 11 shows closed-loop stability at each condition 
(the origin is the critical point). The MIMO Nyquist curves 
are potentially nonconservative measures of stability 
robustness not just for the controller structure shown in 
figure 8, but in general. For rank one loop transfer matrices 
only, as was the case here (see K(s) in figure 8), the Nyquist 
curves directly display tolerances to uniform gain and phase 
errors on all channels or, equivalently, tolerances to gain and 
phase errors at the single strand point. Figure 12 provides 
additional stability robustness information. 

In figure 12 minimum singular value curves &c(o)) and 
£s(o>) are shown for the return difference matrices at the plant 
input and output, respectively. These unstructured singular 
values correspond to the worst combination of independent 
multiplicative errors in the respective input or output channels. 
The magnitude of the MIMO Nyquist curve from figure 1 1 
(plotted as the solid curve in figure 12) corresponds, for this 
rank one controller, to tolerance to uniform (dependent) errors 
and is an upper bound upon Hc(a>) and jLj(co). Finally, a 
minimum singular value curve, Qp(o>), reflecting sensitivity to 
plant additive error, is shown. The singular value assessments 
of stability robustness are potentially quite conservative, since 
the likelihood of encountering the worst combination is not 
addressed. The potentially conservative nature of the singular 
value robustness assessments becomes more pronounced as the 
number of sensor or control channels is increased. 
Nevertheless, the singular values obtained for the MIMO 
controller compared well with those for SISO controllers that 
were tested during this entry, 1 2 » 22 * 23 particularly below 
250 psf. 

All curves for all dynamic pressures show a good balance in 
tolerance to phase lead errors (5.5 < 0 ) < 7.5 Hz) and phase lag 
errors (14 < co < 16.5 Hz). The curves for 290 psf exhibit 
noisy behavior because of a poor signal -to-noise ratio. Only 
small amplitude (±0.3 degree) input excitations could be made 
at 290 psf without exceeding torsional safely load limits. 
Nevertheless, it is apparent that the singular value minimum 



Fig. 1 3 RMS control rates due to turbulence. 
(Symmetric, reoptimized FSS with MLQ 


Command Response 



Fig. 14 Aggressive roll maneuver. 
(11% beyond flutter dynamic pressure) 


levels arc significantly lower than those at the lower dynamic 
pressures. For 290 psf, the lowest singular values are at a 
frequency in close correspondence with the flutter frequency. 
This indicates that somewhat better performance could be 
expected with a higher gain controller. 

Control surface RMS rates during a test run for the free-to- 
roll configuration are shown in figure 13 where they are 
contrasted with analytical predictions. Both actual and 
predicted RMS rates were well within the 50 deg/sec RMS 
design limit. There was qualitative agreement between 
analysis and experiment indicating that the analytically 
assumed gust environment was adequate for assessment of 
potential rate saturation problems. 

Increasingly aggressive maneuvers (in terms of roll 
acceleration commands and dynamic pressures) were tested 
using a multifunction controller (i.e,, both active flutter 
suppression and maneuver load control 18 were simultaneously 
employed). Figure 14 depicts the successful completion of the 
most aggressive maneuver. The maneuver commanded roll 
accelerations of 1000 dcg/sec 2 until a roll rate of 250 deg/sec 
was achieved, held that rate until a 90 degree roll was 
completed (0-to-90 degrees in 0.5 sec), and then decelerated at 
approximately 1000 dcg/sec 2 . The maneuver was performed 
at a dynamic pressure 1 1 percent above the open-loop flutter 
boundary. The aggressive rolling maneuver did not 
appreciably affect the FSS peak control deflection or rate 
requirements. 

A flutter suppression test was conducted with the AFW 
wind-tunnel model in the fixed-in-roll configuration. For this 
configuration, there was considerable disparity between the 
observed open-loop antisymmetric flutter dynamic pressure 
(219 psf) and that predicted with the design model chosen 
(252 psf for the pre-1989 model). Nevertheless, flutter was 
simultaneously suppressed in both symmetric and 
antisymmetric degrees-of- freedom to a dynamic pressure of 
275 psf. This test was terminated at 275 psf because the 
turbulence induced loads were closely approaching preset 
safety limits on torsional loads. The remaining test time was 
too limited to allow generation of data for a full CPE analysis. 
However, stability and RMS control rate usage were 
determined at a number of test points. The RMS control rate 
requirements were well within the design limit being 
comparable to those of figure 13. 

Conclusions 

A multi-input/multi -output flutter suppression system for 
the Active Flexible Wing wind-tunnel model was developed in 
the continuous domain, implemented digitally and tested in the 
Langley Transonic Dynamics Tunnel. Optimization 
techniques were employed to define a controller that exhibited 
significant robustness to flutter frequency prediction errors and 
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to unstructured plant additive errors and unstructured 
multiplicative errors at the plant input and plant output. 

The frequency domain based design approach readily 
permitted utilization of a combination of analytically predicted 
and experimentally estimated plant transfer matrices. An 
overnight reoptimization using combined experimental and 
analytical data was accomplished during the lest. Stability 
robustness improvements resulted as compared with a 
stabilizing design based solely upon analytical predictions. 
The rapid redesign shows the potential, with a much faster 
computer than that used for controller implementation, for use 
of the approach in tuning type adaptive applications. 

The reoptimized control law was tested. Symmetric flutter 
was suppressed to the tunnel limit. With a fixed- in-roll 
configuration, both symmetric and antisymmetric flutter were 
simultaneously suppressed to a dynamic pressure 26 percent 
above the antisymmetric open-loop boundary and 17 percent 
above the symmetric boundary. Turbulence induced loads 
which were approaching a preset torsional loads safety limit 
prevented attainment of the tunnel limit in this case. 
Simultaneous flutter suppression and maneuver load control 
were demonstrated during aggressive rolling maneuvers 
performed at a dynamic pressure 1 1 percent above the open- 
loop flutter boundary. 
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